Calculation of resistance for weak scattering, strong scattering and insulating 

quasi-one dimensional systems 
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A parameter free calculation of the resistivity is applied to liquid metals near the melting point 
ranging from weak to strong scattering limit. The method is based on length dependent resistance 
calculations for quasi-one dimensional systems and was applied on structures with up to 10000 atoms. 
The calculated value for conductance fluctuations is in good agreement with theoretical predictions. 
The resistivities are compared with the Kubo-Greenwood and the extended Ziman formula with 
the same scattering potential and similar structure. The resistance calculation is applicable for 
insulating materials as well, which is demonstrated for crystalline and amorphous silicon. 



I. INTRODUCTION 

Electronic transport in mesoscopic systems is marked 
by interference effects, which can not be understood by 
classical considerations. Effects such as universal conduc- 
tance fluctuations (UCF), Aharonov-Bohm effect, reso- 
nant scattering, etc. have been investigated for years. 
One of the most striking effects is the localization in 
pure one dimensional disordered systems. This led to the 
scaling approach, 1 which demonstrated disorder induced 
transition from metallic to insulating behavior (Anderson 
transition) for three dimensional systems. 

Starting from one dimensional systems^ a scaling ap- 
proach was developed for quasi-one dimensional systems 
employing random scattering matrices^ The length- 
dependence of the scattering system is described by the 
diffusion-type DMPK^ equation, which has a metallic 
solution for lengths shorter then the localization length. 
The resistance increases linear with the sample lengths 
apart from the fluctuations. For large lengths the system 
behaves like a one-dimensional system with a localization 
length dependent on sample size (number of channels) 
and the elastic mean free path. 

The resistance calculations are based on formulas in- 
troduced for one-dimensional systems by Landauer— and 
extended later to quasi-one dimensional systems^. In 
the simplest case the conductance is proportional to the 
total current transmission through the sample, G ~ T, 
which can be understood as "measurement" of the cur- 
rent for a defined voltage drop over the sample. A some- 
what different approach was followed by LenkS using a 
functional of the currents. 

For practical applications the nearly classical scal- 
ing of the resistance was used by Kahnt 10 to calcu- 
late the resistivity for quasi-one dimensional systems of 
muffin-tin scatterer. Subsequently a similar approach 
was used to calculate resistivities within the tight-binding 
formalism^! Both approaches were hampered by large 
deviations from the classical scaling due to the limited 
system size. 

This paper follows the approach of Ref. 0, because it 
is a parameter-free description of real materials. Com- 
pared to previous work the calculation of the scattering 



matrix is done in a combination of plane wave and angu- 
lar momentum representation. This hybrid representa- 
tion preserves the high precision of the angular moment 
representation, which is necessary for resistance calcula- 
tions, but on the other hand avoids the cubic increase in 
the computation time. The computation time increases 
linear with the sample length. Thus up to 10000 atoms 
can be computed with angular momentum channels up 
to I = 2. This makes it possible to extend the resistivity 
calculation on weak scattering systems. 

In section [H] the calculation of the scattering matrix 
and the resistance is outlined. The details of the hybrid 
method are given in appendix [H] 

The sample length dependence of the resistance is ana- 
lyzed for metallic samples in section lTlTl It is shown, that 
there are three different regimes for the length-dependent 
resistance. The theoretical value for the UCF is repro- 
duced. An upper bound for the conductance is found 
for the localized regime. The section is concluded with 
practical aspects of the resistivity calculation, e.g. tem- 
perature dependence. 

The method of resistivity calculation is applied to liq- 
uid metals near the melting point in section llVl Calcula- 
tions for strong scattering 3d transition metals and weak 
scattering liquid sodium are presented. 

Usually the resistivity would be calculated by Kubo- 
Greenwood formulaiSii^ or in case of weak scattering by 
Ziman formula 1 *, therefore the presented method is com- 
pared to these methods using similar input data, e.g. 
structure and scattering potentials. In both cases there is 
a good agreement. Small but systematically larger val- 
ues are found for the transition metals, which may be 
connected to weak scattering corrections. 

For weak scattering much larger samples are necessary. 
It is shown, that the numerical effort can be reduced by 
using an alternative resistance "measurement" employ- 
ing adaptive leads. The resistivity as well as the resis- 
tivity coefficient are calculated near the melting point 
and compared to measurements and calculation with the 
extended Ziman formula 15 using the same structure fac- 
tor and scattering potentials. The influence of multiple 
scattering is discussed. 

The method is finally applied to crystalline and amor- 
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phous silicon. These insulating samples show a different 
length dependence of the resistance, where in contrast 
to localized quasi-one dimensional metallic samples, the 
localization length is characteristic for the insulating ma- 
terial. 



II. CALCULATION OF RESISTANCE AND 
RESISTIVITY 
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A. The scattering matrix of the stack 

The scattering matrix is needed for the calculation of 
the resistance. For most applications the sample length 
dependence of the resistance has to be known. There- 
fore successive calculation of the scattering matrix is em- 
ployed to minimize the computation time. Usually the 
computation is done atom layer by atom layer for a stack 
of layers up to 10000 atoms per sample. 

The scattering of the atoms is approximated by a 
muffin-tin potential calculated self-consistent by a LMTO 
method in local density approximation (LDA)>i& The 
scatterers are described by their energy dependent phases 
shift T)i(E) for the relevant angular momentum I. This 
energy dependence is droped in the following equations, 
because the multiple scattering is calculated for one en- 
ergy E = k 2 . 

Periodic boundary conditions are used in lateral direc- 
tion (perpendicular to the current). Therefore the wave 
field outside is best represented in plane waves 
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exp ik~r 



(1) 



with wave vector 



kf = k\\ +t + e z an? (2) 
and the component in direction of the flowing current 



(fc||+r) 2 , 



(3) 



where a = ± denotes right or left moving waves. Note 
the dependence on the parallel wave vector ku , due to the 
periodic boundary conditions in lateral direction. This 
vector can assume values in a two dimensional Brillouin 
zone due to the lateral boundary conditions. The plane 
wave itself is identified by the two dimensional vector f. 

There are propagating waves and evanescent waves 
with real and imaginary impulse K? perpendicular to 
the stack. Whereas the number of propagating waves 
N « ^—a 2 is finite, the number of evanescent waves is not 
limited. But only a finite number of the evanescent waves 
contribute to the total current, because of the exponen- 
tial decreasing amplitude of these waves. This number 
increases with decreasing distance between neighboring 
layers of atomsiiS 



FIG. 1: Two dimensional outline of the simultaneous use of 
plane waves (pw) and angular momentum (am) representation 
of the wave field. The stack of layers with lateral length a and 
stack length L is increased by an additional sub-stack. 



The scattering of the stacked layers is described by the 
scattering matrix S. The left and right outgoing waves 
are given by the left and right incident plane waves 



r 



S++ S^ 
S"+ S" 



(4) 



where for example 'J; denotes the amplitudes of the in- 
cident and outgoing plane waves on the left of the stack. 
Only the propagating part of the matrix S is used for 
the resistance calculation. These matrix elements are 
the transmission and reflection amplitudes of the stack. 
To distinguish with the scattering matrix in Eq. , the 
transmission and reflection amplitudes will be denoted 
by t and r respectively. Note that only this propagating 
part of the scattering matrix is unitary. 

The successive calculation of the scattering matrix can 
be done easily in this plane wave representation by elim- 
inating the wave field between the main stack and the 
additional layers (Fig. ^). 17 The computation time in- 
creases only linear with the number of layers. This is the 
method of choice, when the number of necessary plane 
waves is small. The plane wave representation is not ap- 
plicable on samples of highly disordered material with 
large lateral length a and high angular momentum chan- 
nels for the scatterers. 

To circumvent this problem the successive calculation 
of the scattered waves was done in angular momentum 
representation in Ref . The propagation between the 
atoms is described by structure constants G l [ L , given in 
appendix 1X1 The scattering is defined by the scattering 
amplitudes 2ikfiL = 1 — e 2lrn , where L = (l,m) denotes 
the angular momentum channel at the i-th atom. The 
wave field in angular momentum basis is given by 



(5) 



where a compact matrix notation is used. The wave field 
is calculated by a successive matrix inversion. 

The plane wave representation for the stack is reached 
by the projection from plane wave to angular momentum 
channels for the incident waves and vice versa for the 
scattered wavesiiS 



S = P hom + ikB [F- 1 - G] *A 



(6) 
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The draw back of the angular momentum representation 
is an increasing number of scattering channels with in- 
creasing numbers of layers. This limits the number of 
atoms to a few hundred for d-scattering. 

As will be shown below, for the calculation of mate- 
rial specific properties, e.g. the resistivity, larger sam- 
ples with more then 1000 layers are necessary to achieve 
a good accuracy compared to Kubo-Greenwood calcula- 
tions. 

Therefore a combination of plane wave and angular 
momentum basis was used to represent the wave field 
in the stacks Whereas far away parts of the stack are 
projected on plane waves, the scattering waves of close 
layers are dealt with in angular momentum basis. This 
limits both the number of plane waves and angular mo- 
mentum channels. The computation time increases only 
linear with the number of layers. The details are given 
in the appendix 151 

B. Calculation of the resistance 

The resistance is usually calculated by a multi-channel 
Biittiker formula 7 

e 2 

G=—T. (7) 

nn 

which connects the conductance G with the total current 
transmission through the stack of layers^ 

The total transmission T is calculated from the trans- 
mission amplitudes Uj of Eq. J3J. Because no scattering 
takes place in the semi-infinite leads to the reservoirs, the 
current flows only through the N propagating channels, 

N N 

r = EEM 2 - ( 8 ) 

i=i j=i 

A shortcoming of Eq. Q is the non- vanishing resistance 
TZq = ^N" 1 , if there are no scattering layers at all. 
The resistance formula Eq. Q is not unique, due to the 
boundary conditions left and right of the stack. The 
ideal leads in connection with the reservoirs define an 
even incident current distribution^* 

To evaluate the influence of the left and right boundary 
conditions on the resistance and more important on the 
resistivity, we compare the ideal leads in case of Eq. Q 
with an ansatz employing adaptive leads;* 

2 N N 
i=l 3=1 

using the elements of the matrix 

T = 2T[1 + R-T] -1 . (10) 

The matrix elements for the transmission Tjy = 

and reflection probabilities By = \rij\ are given by the 

transmission and reflection amplitudes respectively. 




FIG. 2: Comparison of the length-dependent resistance with 
the radial pair correlation function g(r) (- • -) shows the in- 
fluence of the nearest neighbor distance for strong scattering. 

The resistance was calculated with Eq. Q ( — ) and iJS) ( ) 

in T-point approximation at the Fermi energy. The residual 
resistance Hq is omitted for Eq. Q . The lateral length of the 
sample was a = 40 au. Liquid iron near the melting point 
was taken as a typical example for all 3d transition metals. 
The radial pair correlation was taken from Ref. [2oL 

Eq. (JSJ gives a minimal resistance with respect to a 
given functional of the incident currents. 9 In that sense 
it is a good test for the deviations of Eq. J7J). The in- 
verse matrix in Eq. (|10|) characterizes the adaption of 
the currents to the scattering by the stack. 

From Fig. [5] it can be seen, that differences in the re- 
sistance occurred for short stack length. But these differ- 
ences in the resistance remain constant for larger stacks. 
The length dependence is then dominated by the trans- 
mission. Nevertheless the resistance in Fig. [5] is always 
smaller for adaptive leads then for ideal leads. 

III. RESISTANCE FOR DISORDERED 
METALLIC SAMPLES 

A. Dependence on the sample length 

For a disordered metallic sample the length depen- 
dence of the resistance can be divided in three different 
regions. 

For a short length the stack is not material specific. 
This can be seen from the difference of the measured 
resistance, i.e. Eq. Q and ©. This not material specific 
region is given by the characteristic length scales of the 
material, e.g. the correlation length and the elastic mean 
free path. 

From Fig. [21 most strikingly visible is the connec- 
tion to the interatomic distance. This is caused by the 
strong scattering between nearest neighbors in 3d tran- 
sition metals. This strong scattering enables an effective 
current transmission through the sample, if atoms are 
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FIG. 3: Variance of the conductance for a = 40 au and a — 20 
au (inlet) plotted over the sample length L. The theoretical 
value is 2/15^ 50 samples as shown in Fig. |5|were used for 
a = 40 au and all the samples for the inlet are shown in Fig.0 



stacked one behind the other. Therefore the slope of the 
resistance is usually smaller beyond this initial region. 

The second region is material specific. The resistivity 
can be calculated from the slope of the resistance over 
the sample length. This is done by linear regression and 
can be seen as a fit to the classical relation 



TZ = K b + p 



L 
A 



(11) 



for a metallic wire of length L and cross section A = a 2 . 
The initial region discussed above is accounted for by a 
length independent resistance TZf,. 

As can be seen from Fig. there is no notable change 
in the difference between the resistance calculated with 
Eq. and Eq. © in the metallic region. Therefore 
the linear regression with Eq. (jnj gives a "measure- 
ment" independent resistivity. Length dependent devi- 
ations for different "measurements" of the resistance in- 
dicate a non-material specific behavior (see sec. IIV Bfl . 

The deviation from the classical behavior of Eq. is 

2 

of order ^ for the conductance and known as Universal 

Tin _ _ _ 

Conductance Fluctuations (UCF). In Fig.|3|the variance 
of the conductance is plotted over the sample length L 
for two different lateral length a. The smaller sam ples 
were calculated by the original method used in Ref. [HJ. 
The hybrid representation of appendix iBl has to be used 
for the doubled lateral length. 

For sample length larger then the lateral length, i.e. 
L > a, there is a good agreement with the theoretical 
value 2/15 for UCF in quasi-one dimensional systemsi^Si 
This agreement shows that the variations in the resis- 
tance curves stem from the coherent scattering and not 
from variations of the underlying material. This is im- 
portant for calculations of material specific resistivities. 

From the practical point of view, the conductance fluc- 
tuations result in variations of the resistance slope, most 
notably for small later length. The resistivities for the 
smaller samples of Fig. [3] lie between 65 and 145 /ificm, 



which makes a calculation of a single sample to unreliable 
for comparison with experiment and Kubo-Greenwood 
method. 

By doubling the lateral length the variation of the re- 
sistivity is reduced below ±10 /^f2cm. Thus the material 
is much more defined and the resistivity can be calcu- 
lated from one sample. An averaging over many samples 
is not necessary. 

Anyway, if the length L increases the conductance be- 
comes of the same order as the fluctuations. In Fig. 0] 
this is shown for many samples with small lateral length. 

The average of the conductance decreases exponen- 
tially. Different averaging schemes have been proposed 
and tested in Ref. I 111 The quasi-one dimensional samples 
become localized due to the disorder just as disordered 
one dimensional systems for large L . 

From Fig. |4| it can be seen, that G = ^- is a lower 
boundary for the localized samples. If the total trans- 
mission of a specific sample falls below T = 1, then it 
does not increase beyond that limit for any larger length. 
The wave field is characterized by the localization length 
l\ oc « Nl e \j* where Z e i is the elastic mean free path. The 
resistance is not solely depending on the material like in 
the second regime, but also on the lateral constrains and 
is therefore not material-specific. 

In this localized regime the current is dominated by 
transmission through isolated quasi-bound states in con- 
trast to the metallic regime, where many states con- 
tribute to the current. If the resonance energy coincides 
with the calculated energy and the quasi-bound state lies 
in the middle of the stack, one gets a total transmission 
T = 1. In any other case the transmission is lower. 

It remains to note, that one should be careful with av- 
eraging over different probes. The inherent fluctuation 
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FIG. 4: The inverse average of the conductance (G)~ calcu- 
lated with Eq. J7) is plotted for long sample length L. The 
resistance of all 24 samples are plotted with gray lines. The 
boundary between metallic and localized region (G — ■%■) is 



marked by (- 



The transition region, when the first and 



the last sample falls below this threshold (G = ^r), is out- 
lined by vertical lines (comp. Fig.|^J. The liquid iron samples 
had a lateral length of a ~ 20 au (N = 20). 
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FIG. 5: Part of Fig. |I]is shown with linear scaling of the re- 
sistance. The start of the transition to localization is marked 
by the vertical line. The linear extrapolation of Eq. ill 1 ^ is 

plotted as ( ). The fit was done for the inverse averaged 

conductance in the interval (20 . . . 50) au. 

of the resistance can be reduced by averaging, but only 
in so far as the resistivities of the different probes are the 
same. This problem becomes apparent in the broad tran- 
sition region to localization marked in Fig. Averaging 
of the resistance gives the impression of a non-linear de- 
pendence of the resistance on L. This is best met by an 
average over the conductance (see Fig. because this 
corresponds to a parallel set-up of the probes. Compared 
to the averaging of the resistance, samples with higher 
conductance are weighted stronger. One should note 
however, that a sample average for transmission above 
and below the threshold T = 1 does mean averaging over 
different 'materials', i.e. metallic and insulating samples. 
As can be seen in Fig. this leads to deviations from 
the linear resistance scaling of Eq. . 

This is the reason, why the resistivity calculations 
shown later rely on one larger sample, rather then many 
small ones. Even though the hybrid representation out- 
lined above can be used for calculation of large stack 
length, the main advantage is to increase the lateral 
length and thereby to reduce the difference between sam- 
ples in terms of the resistivity. Just as for standard cal- 
culation of the resistivity by Kubo-Greenwood formula, 
one has to compute only one sample. 

B. Calculation of the resistivity 

Due to the periodic boundary conditions in the lateral 
direction the conductance of Eq. 10 and Eq. © depends 
on the fey-point in the resulting two dimensional Brillouin 
zone. Usually the calculations of the resistance is done 
in T-point approximation, that is only one point in the 
middle of the two dimensional Brillouin zone is used to 
calculate the current through the sample. 

The error made by the single fen -point approximation 
is illustrated in Fig. for a small lateral length and in 



Fig. for the doubled lateral length. The differences in 
the resistivity decrease by increasing the lateral length 
a. For small a the resistivity values differ in the same 
order of magnitude as for different samples in the T-point 
approximation. 

In general the conductance has to be averaged over all 
/c||-pointsiS 

G = ^ |d 2 fc||G(fc||) (12) 

BZ 

This fey -dependence is in so far important, that there are 
cases, when large lateral length a are not attainable and 
the accuracy is not sufficient. An example is a small 
approximants of a quasi crystal,— where the structure 
is well defined. Apparently an averaging over samples is 
not possible. In such a case one has to use Eq. (fT2"|) . 

It should be noted however, that contrary to different 
samples, the dependence on fey in Fig. is far from ran- 
dom. Eventually the resistivity is dominated by small 
regions with small resistivity. Thus one has to calculate 
many different -points to get an accurate result, which 
makes this sampling more time consuming then a calcu- 
lation of one large sample in T-point approximation. 

Both examples shown here have the maximal value of 
the resistance near the far edges of the BZ. Whereas the 
T-point gives a fair approximation of the averaged resis- 
tivity based on Eq. I|12l) . The k\\ -dependence is smoother 
for the larger lateral length. Hence a smaller number of 
sampling points is sufficient to approximate the integral 
in Eq. lfT2}. 

For the calculations corresponding to Fig. the resis- 
tivity is 83 /xficm, compared to 92 /ificm for the larger 
sample of Fig. Note however, that other samples with 
a = 20 au showed up to 15 /ii7cm higher resistivities. 

To conclude this section we note, that for temperatures 
different from zero the conductance has to be averaged 
over energies different from the Fermi energy. Due to 
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FIG. 6: Resistivity of liquid iron calculated for single fcy- 
points of the two dimensional BZ based on Eq. 0. The 
contour lines are plotted every 10 /if2cm starting at 70 ^iilcm. 
The lateral length of the sample is a = 20 au. 
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FIG. 7: Resistivity of liquid iron calculated for single Ap- 
points as shown in Fig.|S|but with a lateral length a = 40 au. 
The contour line is plotted for 90 /jficm 



the phase shift of the scatterers and the propagation of 
the waves in the stack, the conductance depends on the 
energy. Different energies have to be taken into account 
especially for higher temperatures. This results in an 
average of Eq. Q weighted by the Fermi distribution. 23 
Whether calculations for different energies are necessary 
depends on the change of the current transmission in an 
energy interval of order k^T around the Fermi energy. 
For the liquid transition metals the conductance does not 
change much on the energy scale k^T and the average can 
be approximated by 



and experimental measurements from Ref. l25i The re- 
sistance calculations are done in T-point approximation 
by Eq. . Because of the high melting temperatures for 
these metals, five different energies were used to approxi- 
mate Eq. (|13fl . Notable deviations of resistivity from the 
value at the Fermi energy were observed for liquid tita- 
nium and nickel. This is caused by the strong change 
in the DOS around the Fermi energy for nearly filled or 
empty d band. 

All stacks had a lateral length of 40 au and the linear 
regression was done between 10 au and 100 au. In favor of 
a much faster calculation, a full projection on plane waves 
was done only after every 80th layer. The hybrid basis 
consisted of 9 angular momentum channels for each of the 
last 80 attached layers and about 320 plane waves. Both 
values were tested by increasing the number of channels 
and by the accuracy of the current conservation. 

The resistivities calculated by linear regression of the 
multi-channel Biittiker formula Eq. (JJJ are compared 
with Kubo-Greenwood resistivities in Fig. [H] For com- 
pleteness the experimental values are shown as well. The 
differences between numerical and experimental results 
are discussed in Ref. |3 Essentially this is caused by 
the lack of spin dependence in the calculations. Thus 
the differences to experiments are larger for half filled d 
band. 

To demonstrate, that the agreement with Kubo- 
Greenwood resistivities is not limited to dominating d- 
scattering, the resistivity of liquid copper is added in 



G(T) 



dE 



dE 



G{E) 



(13) 



As an example the temperature dependent resistivity and 
the resistivity coefficient is calculated for liquid sodium 
m 

section EH 



IV. APPLICATION ON DIFFERENT 
MATERIAL CLASSES 

A. Liquid 3d transition metals 

In the following section results of resistivity calcula- 
tions are compared with standard methods, i.e. Kubo- 
Greenwood and Ziman formula. To ensure that varia- 
tions in the results stem from the methods not from the 
input data, the same self-consistent LMTO program was 
used for the calculation of the scattering potential^ Be- 
cause the super cells are usually much smaller then the 
stacks, the potential has to be averaged over all atoms of 
the same type. 

The structures were generated by Reverse Monte Carlo 
method based on measured structure factors for the melt- 
ing points The structure factor was matched closely by 
the RMC method aside from liquid titanium. 

In Fig. [S] the resistivities of liquid 3d transition met- 
als are compared with Kubo calculations form Ref. |^ 



200 



Ti V Cr Mn Fe Co Ni Cu 




FIG. 8: Resistivity for liquid 3d transition metals and cop- 
per near the melting point. The calculations based on 
Eq. Q (o) are compared with Kubo-Greenwood results (□) 
and measurements (•)■" Values from left to right in /ificm: 
159.7; 112.7; 97.9; 99.3; 90.8; 98.4; 91.5; 28.3. 
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Fig. [5] Both numerical calculations are in good agree- 
ment with experimental values^ 

Agreement between the resistivity based on Eq. 111(1 
and the Kubo-Greenwood resistivity is remarkably good 
for all liquid 3 d transition metals. This shows that both 
methods are equally suited for resistivity calculations of 
strong scattering disordered materials. The only notable 
deviation is the up to 10 fiflcm larger resistivity for the 
results by multi-channel Biittiker formula. Likewise there 
is no notable difference between the resistivities based on 
Eq. J7J) and Eq. 0. Both calculations use similar struc- 
ture and scattering potentials, periodic boundary condi- 
tions in lateral direction and are based on a linear re- 
sponse ansatz. Likewise the incident current distribution 
has no effect on the resistivity. The only two differences, 
which may have an effect on the resistivity, are the dif- 
ferent lateral dimension and the attenuation used for the 
Kubo-Greenwood calculation in Ref. I24I 

A similar difference for the resistivity occurred for the 
smaller samples used in section IIII Al in case of liquid 
iron. It should be noted, that the size of the supercells 
for Kubo-Greenwood calculations was about the same as 
the lateral length a ss 20 au of the smaller samples. 

These differences in the resistivities resemble the quan- 
tum mechanical corrections due to multiple scattering in 
three dimensional finite systems^ 

Hn 3 \l e i a) ' ^ ^ 

where a is the size of a cubic sample and l e \ is the elastic 
mean free path. 

For the example of liquid iron the Joffe-Regel limit is 
reached implying w 5 au. 24 This gives a weak localiza- 
tion correction for the resistivity of about 8 /zficm, which 
is of the same order compared to the 7 /iQcm higher resis- 
tivity for liquid iron in Fig. EI The increase in resistivity 
for the doubled lateral length a is 1 /iSlcm according to 
Eq. (|14|l and is therefore much too small to account for 
the differences in Fig. [S] 

Hence the main reason for the discrepancies seems to 
be the attenuation of the scattered waves in Ref. l24l 
which may have an effect on the weak localization cor- 
rection. One should note however, that the discussed 
discrepancies are in the range of the error margin espe- 
cially for small lateral length. 

B. Application for weak scattering 

Previous resistivity calculations based on Eq. Q were 
limited to fairly strong scattering, because of the limited 
sample size»iS This constraint is relaxed by using the hy- 
brid representation. 

Different resistance calculations for liquid sodium near 
the melting point are shown in Fig. [§J The transient 
region is much larger then for strong scattering systems 
(comp. Fig. |2J). The lateral length of 70 au was chosen 
according to the initial transient region for the resistivity 




L/[au] 

FIG. 9: Resistance over sample length L for two different lat- 
eral length (a = 50, 70 au). The sample with the smaller lat- 
eral length of 50 au is marked by ♦. Two different resistance 
formulas are used, which correspond to ideal leads Eq. Q and 
adaptive leads Eq. © leads respectively. The latter resistance 
values, i.e. Eq. @, are plotted as broken lines for both lateral 
lengths. Calculations were done in F-point approximation. 



calculations. The results was tested by samples with a = 
50, 100 au. 

Whereas for strong scattering there are only small dif- 
ferences between the resistance formula Eq. Q using 
ideal leads, and Eq. ©, the deviations become practi- 
cally important for weak scattering. For adaptive leads, 
i.e. Eq. ©, the transient region is smaller, because the 
incident current distribution is already adapted to the 
scattering of the sample. 

This adaption has to take place inside the sample for 
Eq. (7J. Due to the weak scattering, long samples are 
needed to reach a current distribution independent from 
the even distribution of the incident currents. 

The slope of Eq. {7|) in the transient regime is larger 
then in the metallic regime. For instance a cubic cell 
with length a = 30 au would give a resistivity of about 
45 /xficm. 

Although not comparable, Kubo-Greenwood calcula- 
tions with similar small cells give far to high resistivi- 
ties near the melting point ^ whereas calculations based 
on Ziman's formula give good results near the melting 
point 

Increasing the lateral length increases the conductance 
per cross section based on Eq. J7J), whereas values for 
Eq. 10 are nearly independent of the lateral length even 
in the transient region (see Fig. . Both equations have 
similar transient regions for a lateral length of a = 70 au. 
Kubo-Greenwood calculation with such a cell dimension 
should give similar results for the resistivity. 

All the shown resistivities are calculated using Eq. , 
because of the smaller dependence on the lateral length 
a (see Fig.0. 

The resistivity calculated for single k\\ -points differ by 
about 2 /ificm for a — 70 au (3 /^f2cm for a — 50 au). 
Sampling was done with 16 different -points to approx- 
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TABLE I: Calculated values of the resistivity near the melt- 
ing point for liquid sodium are compared to Kubo-Greenwood 
formula (KB) and Topp-Hopfield (TH) pseudopotential, 27 
Ziman formula with Shaw's pseudopotential (TW)[— and 
measurementsi— The present calculations are based on the 
same self-consistent calculated LMTO potential. The ex- 
tended Ziman formula— is compared with linear regression 
results for Eq. HI IB based on Eq. © (Lenk). The higher an- 
gular momentum scattering was disregarded for the last two 
calculations (s). 



T(K) 


p (pQcm) 


method 


potential 


400 


~ 25 


KG 


TH 


373 


8.6 


Ziman 


TW 


378 


9.5 


ext. Ziman 


LMTO 


378 


9.2 


Lenk 


LMTO 


378 


10.1 


ext. Ziman 


LMTO (s) 


378 


9.1 


Lenk 


LMTO (s) 


373 


9.44 


Exp. 





imate Eq. I|12l) . The difference to samples with a = 50 
au is reduced below 0.4 fiQcm. Such a high accuracy is 
necessary for the calculation of the resistivity coefficient 
near the melting point 

^ 1 p(T) - p(T m ) 
a ~ (rr \ — ™ — ™ • (15) 

It can be seen from table that calculations based 
on Eq. (PJ are close to the resistivities based on Ziman 
formulai22i22i2i To underline this, values for the resistiv- 
ity are included using the extended Ziman formula of 
Evans— i Because the same phase shift and structure fac- 
tor is used, the deviations are caused by multiple scat- 
tering. 

As for strong scattering the sample structures were de- 
termined by RMC method based on measured structure 
factors^ The phase shifts (up to I = 2) were calculated 
self-consistently by a LMTO method^ Half the mini- 
mal distance between atoms was chosen for the muffin 
tin radius. 

Naturally s-scattering is dominating the resistivity for 
liquid sodium, but nevertheless the inclusion of p- and d- 
scattering yields to a notably reduction of the resistivity. 
This effect is small near the melting point and is nearly 
canceled due to multiple scattering (see Tab.HJ, but be- 
comes more important for higher temperatures (comp. 
Tab.HU). 

The calculated resistivities based on the extended Zi- 
man formula as well as multiple scattering are in good 
agreement compared to former calculations based on the 
Ziman formula,— and in even better agreement to the 
measured resistivities near the melting point ^2i2i But so 
far the Ziman formula failed to describe the temperature 
dependence— i 

As can be seen from Tab.lTTl the deviations from mea- 
sured resistivities increase with increasing temperature, 
because the calculated resistivities remain comparatively 
unchanged. Nevertheless, there is the good agreement 



TABLE II: Values of the resistivity for higher temperatures 
compared to measured resistivities from Ref. l30t Abbrevia- 
tions are the same as for Tab. [I] 



T (K) 


p (pQcm 


) method potential 


473 


11.5 


ext. Ziman LMTO 


473 


10.5 


Lenk LMTO 


473 


13.95 


ext. Ziman LMTO (s) 


473 


11.8 


Lenk LMTO (s) 


823 


11.7 


ext. Ziman LMTO 


823 


11.4 


Lenk LMTO 


423 


11.1 


Exp. 


473 


12.9 


Exp. 


823 


28.56 


Exp. 


TABLE III: Values of the resistivity coefficient near the melt- 


ing point. 








Exp. 


Ziman ext. Ziman Lenk 


a ■ K/IO* 


3.5 


2.9 2.2 (4.0) 1.5 (3.1) 



between the multiple scattering results based on Eq. (jl 1|) 
and 10 and the extended Ziman formula. 

The resistivity coefficients are summarized in Tab. Mil 
The resistivity coefficients are lower compared to the ex- 
tended Ziman formula, if multiple scattering is included. 
This is the same behavior as for the resistivity. Likewise, 
the resistivity coefficients are increased, if higher angular 
momentum channels are excluded. 

The results are comparable to the value based on mea- 
surements for the reduced number of channels. Possible 
reasons are relative inaccurate phase shifts for higher an- 
gular momentum channels, because these phase shifts are 
quite small. A hint is the strong dependence of the calcu- 
lated phase shifts for higher angular momentum channels 
on the muffin tin radius. 

Anyway the agreement may be accidental, because 
both the extended Ziman formula and the presented 
model fail to describe the temperature dependence of the 
resistivity over a larger temperature range (see Tab. ILTI). 



C. Resistance for insulating materials 

So far it was shown, that there exists a region, where 
the resistance scales linear with the sample length (apart 
from UCF). For large sample length, there is a transition 
to localization due to disorder. This transition is caused 
by the lateral confinement. Thus the localization length 
^loc ~ Nl c \ depends on the number of propagating modes 
N and the elastic mean free path l e \& 

In comparison the resistance is shown for samples of 
crystalline and amorphous silicon in Fig. 1101 The resis- 
tance increases exponentially in both cases, in contrast 
to metallic materials (see Fig. QJ, where an exponential 
increase is only observed for the sample average. 

The exponential scaling of the resistance is indepen- 



FIG. 10: Resistance of amorphous ( — ) and crystalline silicon 
samples (...) in T-point approximation in the middle of the 
gap (see Fig. II IB . The lateral length of the samples is about 
20 au for a-Si and c-Si. An additional sample with a ~ 10 
au was calculated for c-Si, but it is indistinguishable from the 
sample with a m 20 au in this representation. 



dent from the lateral length of the sample. The localiza- 
tion length is characteristic for the insulating material of 
the sample in contrast to metallic samples, where the lo- 
calization length increases with increasing cross section. 

Modelled structures show usually a distinctive metal- 
lic regime, because the energy of the electrons is larger 
then the scattering potentials. Therefore the exponen- 
tially decreasing current transmission has to result from 
a quite effective destructive interference of the scattering 
waves. Essentially these interferences result in a reduced 
density of states (DOS) near the Fermi energy, which is 
connected to the diagonal parts of the inverse matrix in 
Eq. ©pUii The resulting gaps for crystalline and amor- 
phous silicon are shown in Fig. ^2 Note however, that 
there is no one to one correspondence between electronic 
defects in Fig. 1111 and resonances of the transmission in 
Fig. El because of the different employed boundary con- 
ditions and calculation methods. 




-1.0 -0.5 

(E-E F )/[Ry] 

FIG. 11: Electronic density of states for a-Si ( — ) and c-Si 
( ) calculated with a self-consistent LMTO method. 



FIG. 12: Energy-dependent conductance of the a-Si ( — ) and 

c-Si ( ) samples shown in Fig. 1101 The length L is about 

160 au for both samples. 



Both DOS calculations in Fig. ^2 were done with a 
standard self-consistent LMTO-method in atomic sphere 
approximation. Only s- and p-electrons were used for the 
amorphous silicon calculation, because the cell included 
512 silicon atoms (compared to 8 for the c-Si). Vacancies 
had to be included in the calculations to get better results 
for the DOS, since both silicon modifications are open 
structures. About half the number of silicon atoms was 
sufficient to get a reasonably filled structures: 39 

The amorphous silicon networks were generated by 
a RMC method with subsequent relaxation by a MD 
method. The first step was introduced to ensure an amor- 
phous network structure with four fold coordination. The 
network configuration was fitted to the radial pair distri- 
bution based on measured structure factor^ The latter 
step was necessary to get the gap at the Fermi energy 
and is based on multiple scattering. 33 

The combination of both methods for generating amor- 
phous silicon structures will be described elsewhere. The 
DOS is comparable to ab initio results^ The stabil- 
ity of the final structures was tested by an ab initio 
method (ABINIT) for a smaller sample with 64 silicon 
atoms' The energy dependent transmission through 
such an amorphous silicon network resembles the trans- 
mission through the crystalline counterpart at least in the 
middle of the gap (see Fig. ITUland (|El|) ). Electronic de- 
fects in amorphous silicon are visible in the conductance 
especially for the tail region of the DOS (see Fig. I12J) . 

This part in the energy dependent conductance plot 
increases for larger sample length. A conductance plot 
over the sample length in this energy region corresponds 
to the localized part of Fig. 0] where the conductance 

2 

is lower then Only where the total transmission is 
lower then one, the energy dependent conductance can 
be represented by a one dimensional model of coherent 
scattering. 
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V. CONCLUSIONS 

Starting from scattered-wave calculations of meso- 
scopic quasi-one dimensional systems, the length depen- 
dence of the resistance can be divided in three different 
regimes. For short length, the calculated resistance de- 
pends on the employed resistance formula. For larger 
length the increment of the resistance becomes indepen- 
dent of the incident currents and is characteristic for the 
material. 

This second metallic regime can be used to calculate 
the resistivity for a wide range of disordered materials, 
because the resistance shows a nearly classical depen- 
dence on the sample size. It was shown, that the re- 
sults are in accordance with standard methods like Kubo- 
Greenwood and Ziman formula. Likewise there was a 
reasonable agreement with measurements. A large num- 
ber of atoms can be calculated with the presented hybrid 
representation of the wave field to get a high accuracy 
for the calculated resistivity. The precision is further in- 
creased by sampling over different fey -points. 

The resistance shows large fluctuations for large sam- 
ple length due to disorder induced localization of the scat- 
tered waves. The sample averaged resistance increases 
exponentially. The total transmission through the stack 
is lower then one. The length dependence resembles a 
disordered one dimensional system. 

For another coherence effect, the Universal Conduc- 
tion fluctuations (UCF), it is shown, that the theoretical 
predicted value 2/15 for a quasi-one dimensional system 
is met quantitatively. This can be used as a test, how 
well defined the material properties are. 

For the example of strong scattering liquid transition 
metals a somewhat larger resistivity was observed, com- 
pared to Kubo-Greenwood calculations. This was inter- 
preted by weak localization corrections known from scal- 
ing in three dimensional systems. This underlines the 
similarities between the scattering in (finite) three dimen- 
sional systems and in the metallic regime of a quasi-one 
dimensional system. 

As a second example the resistivity of weak scatter- 
ing liquid sodium near the melting point was calculated 
and compared to extended Ziman formula using the same 
structural and scattering properties. The agreement was 
good, but multiple scattering corrections seem to be im- 
portant especially for the resistivity coefficient. Both the 
resistivity and the resistivity coefficient are lowered by 
multiple scattering. 

Whereas the calculated resistivity near the melting 
point is close to measured values and even the resistiv- 
ity coefficient may be explained, measured resistivities 
for higher temperatures are not matched by the present 
calculations. 

As a final application of the method, the transmission 
through insulating crystalline and amorphous silicon was 
calculated. The atomic models show an exponential de- 
crease of the transmission like in case of a simple one 
dimensional potential barrier. 



The energy dependent transmission through the barri- 
ers was compared. Differences between amorphous and 
crystalline samples occur at the band edges, where the 
disorder leads to quasi-bound states. This region corre- 
sponds to the localized regime observed for length depen- 
dent resistance of metallic samples. 

APPENDIX A: PROPAGATORS 

1. angular momentum basis 

The propagation of the waves in angular momentum 
representation between muffin tins in different layers is 
described by the interplanar structure constants 

G'lv = £C(£|£|I/)G* = J2 C ( L '\ L \ L ) G oL > ( A1 ) 

L L 

with 

G Lo = E ikh ^ {fi " ■ (A2) 

The propagation between muffin tins inside a layer is 
described by 

G' w = £ C(L\L\L')G' Lo = C(L'\L\L)G' oL (A3) 

L L 

with 

G' L0 = E ikh<»(rj-i*p)e**<&^, (A4) 

and C{L'\L\L) the Gaunt coefficients. A direct calcu- 
lation of Eq. (|A1|) and (|A3fl . as well as a computation 
in reciprocal space is not possible, therefore an Ewald 
summation is used^& 



2. Transformation between plane waves and 
angular momentum representation 

Transformation from incident plane waves in angular 
momentum basis gives the matrix elements 

(A). Lfa = -^^(-l)- e ^(^--«)y_ i (%) . (A5) 
v K f — 

The back transformation of the scattered waves is given 
by 

(B) r - QjL = 1 -^ I V4^e ik ^-^Y L (kl) . (A6) 

The propagation of electrons without scattering is best 
described in plane wave representation 

/phom\^ r ,X^,p&"(fr-i f i) ( \J\ 

where r± — ry /; are the reference points right and left of 
the stack. The derivation can be found in Ref. |37J. 
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APPENDIX B: MULTIPLE SCATTERING IN 
HYBRID REPRESENTATION 

The projection of one part of the stack on plane waves 
basis 

\ L t>\ nc ) = Ai |*\pj nc ) , = Bi \ L $l c ) (Bl) 

is the starting point of hybrid method. The representa- 
tion is distinguished by the index p for plane wave and 
L for angular momentum basis. For convenience the re- 



maining part of the wave field is divided in the notation. 
The scattered wave for the initial step is given by 

/ \ / Tg T?£ T\l \ I \ 

= T % T LL T LL L WT c (B2) 

V L *r / V t& / V L n nc J 

The additional substacks are attached by increasing the 
angular momentum basis 



Lysc 



+ 



/ rppp rppL 

I - L ll - L 12 



rpLp rr\LL rpLL 

A 21 ± 22 s - r > 



3 

V o 

/ 



-32 









\ rpLLrjLp 
\ J- 44 f 41 



-23 



TLp rr\LL rpLL 
31 - 1 - 39 J- 3 



-33 







rrLLr^LL 

44 -^42 



^ 




/ P\J7*™ C \ 














L^mc 


x 44 / 











rpLLryLL 
*-4A *43 



± 11*14 
rpLpr,pL 
*-21 r 14 
rpLp-ppL 
A 31 r 14 



' x 12 *24 ' 

TLLrtLL 
22 *24 

-■-32 *34 





rpPi-pLL 

' x 13 *34 

" - L 23 ^34 
rr^LLnLL 

- X oq JT o/i 



\ 


( \ 













(B3) 



where T 44 describes the scattering by the isolated sub- 
stack and 

■ppL _ fpfiomp ] — -pLp _ [* pfcoml + 

P m4 = , P 4r ^ = Gf m , to = 2, 3, (B4) 

are the propagators between the old main stack and the 
sub stack (s. appendix El . The index a — ± indicates 



the direction of propagation. 

In the third step Eq. I|B3|) is solved, which results in a 
hybrid-matrix R, 

|* sc ) =R|* inc ) . (B5) 

In the last step, parts of the resulting hybrid-matrix are 
projected on the plane wave basis 









-( 







Rff + R#A ? + B 2 (R£f + R^A 2 ) R^g + B 2 R 23 L + B 2 R£ 4 L \ / 



R| + RgA 2 Rg Rg I ' « L5(ii 

R 4 f + R 42 L A 2 R-43 L ^44 / V ^4 



and the right reference point is moved to the right of 
the stack. The hybrid representation of the increased 
stack has the same structure and dimension as the initial 
stack (see Eq. (|B2|l '). Projection of the remaining parts 
on (propagating) plane waves gives Eq. (jSJ. 

Repetitive applying of these four steps gives the wave 
field of the successive increased stack without increase of 



the numerical cost per step. The accuracy were tested 
by current conservation. Stacks with up to 10000 d- 
scatterers were calculated. The in general increasing 
error in the current conservation is not dominated by 
the above inversion but by the Ewald summation for the 
structure constants of appendix lAl 
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The lactor 2 for the electron spin is already included. 
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39 An equal number of silicon atoms and vacancies were used 
for c-Si calculations. 



